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In  factorial  experiments  the  generation  of  treatment  combination  is  the 
beginning  of  the  statistical  analysis.  Plans  presented  in  standard  texts 
are  not  necessarily  suitable  in  most  cases.  A  systematic  procedure  for 
generating  treatment  combination  has  been  presented  by  [1  ]  and  latter 
was  extended  by  [2],  In  this  research  we  show  an  alternative  method 
using  Pascal  triangle  and  lexicographical  order  to  obtain  number  of 
factorial  effects  and  treatment  combinations  respectively.  The  new 
procedures  save  space  and  time. 


1.  Introduction 

Design  of  experiments  has  played  a  fundamental  role  in  the  statistical  curriculum,  practice,  and 
research  ever  since  R.A.  Fisher  founded  the  modem  discipline  [3].  It  has  been  successfully  applied 
in  many  fields  of  scientific  investigation.  These  include  agriculture,  medicine,  and  behavioral 
research  as  well  as  chemical,  manufacturing  and  high-technical  industries.  Concepts  like 
randomization,  confounding,  and  aliasing,  which  originated  in  factorial  experiments,  have  found 
applications  beyond  their  initial  motivation.  An  important  application  of  statistical  methods  to 
industrial  research  is  in  the  factorial  experiments  in  connection  with  the  improvement  of 
manufacturing  processes.  This  applies  to  chemical  and  biological  processes,  formulation  of 
pharmaceutical  preparations  and,  in  fact,  to  most  types  of  industrial  research.  For  example,  in 
chemical  research  it  is  frequently  required  to  determine  the  effect  of  certain  changes  in  reaction 
conditions  or  methods  of  manufacture  on  the  yield  and  quality  of  chemical  products.  Typical 
reaction  conditions  which  might  be  varied  are  temperature  of  reaction,  time  of  reaction,  rates  and 
methods  of  agitation,  concentration  and  amounts  of  reactants,  different  catalysts,  etc.  Similar 
considerations  arise  in  chemicals  manufactured  by  fermentation  processes,  e.g.  penicillin, 
streptomycin,  industrial  alcohol,  lactic  acid,  etc.,  for  which  we  may  wish  to  examine  the  effect  of 
changes  in  the  conditions  of  fermentation  on  the  yield  and  quality  of  the  products.  The  objective 
of  all  such  investigations  is  usually  to  improve  the  yield  or  the  quality  of  the  product,  or  to 
produce  the  product  more  economically. 

The  factorial  experiment  in  which  each  of  the  possible  combinations  of  the  levels  of  several 
factors  are  allocated  to  one  or  more  experimental  units  is  an  experimental  technique  [4].  An 
alternative  formulation  for  such  treatment  combination  was  also  given  using  Kronecker  product 
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[5].  Other  standard  text  books  that  also  reported  the  method  of  Yate  order  include  [6,  7,  8,  9,  10], 
Hand  12]. 

Factorial  design  is  a  strategy  in  which  factors  are  simultaneously  varied,  instead  of  one  at  time 
experiment  [13].  The  common  factorial  design  is  2n  which  is  applicable  when  many  factors  are  to 
be  investigated,  with  the  aim  of  finding  out  which  factorial  effects  and  interaction  between  factors 
are  the  most  influential  on  the  response  of  the  experiment  where  n  is  the  number  of  factors. 
Factorial  experiments  have  come  to  play  a  prominent  role  in  such  investigations  because  of  their 
advantages  in  economy,  precision  and  accuracy.  Besides  having  applicability  to  the  screening  of 
factors  and  the  study  of  the  ranges  of  factors,  i.e.,  the  exploratory  stages  of  scientific 
investigation,  factorial  designs  are  of  use  in  studying  detail  qualitative  factors  such  as  varietal 
differences. 

Factorial  designs  have  been  used  for  many  years  in  agricultural  and  biological  research,  where  the 
practice  has  been  to  produce  a  large  number  of  observations  relative  to  the  number  of  factors 
studied  (e.g.,  by  replication). 

An  early  stage  in  the  construction  of  two-level  factorial  experiments  is  the  generating  of 
treatments  combinations  and  besides  the  Yates  standard  order  as  proposed  by  [14],  most 
researchers  and  standard  texts  give  no  guidance  as  to  the  generation  of  these  treatment 
combinations  other  than  by  trial  and  error,  hence  there  is  need  to  come  up  with  a  systematic 
procedure  to  carter  for  this  problem. 

1.1  Symmetric  Factorial  Designs 

Although  the  2n  factorial  represents  only  a  special  case  of  the  general  Sn  factorial  experiment, 
that  is,  n  factors  with  S  levels  each,  it  deserves  special  consideration  because  of  its  practical 
importance  and  its  special  algebraic  and  combinatorial  representation.  The  list  of  treatment 
combination  is  called  the  design  matrix  and  is  denoted  by  D.  For  a  2n  factorial,  the  design  matrix 
contains  n  columns  and  N  =  2n  rows.  There  is  a  column  for  each  of  the  n  variables,  and  each  row 
gives  a  treatment  combinations.  Normally  the  runs  are  listed  in  standard  order  due  to  [1],  An 
experiment  in  which  each  of  the  2n  possible  combinations  of  factor  levels  occurs  as  a  treatment  is 
called  a  complete  2n  design.  The  two  levels  for  each  factor  are  generally  called  low  and  high. 
These  terms  have  clear  meanings  if  the  factors  are  quantitative,  but  they  are  often  used  as  labels 
even  when  the  factors  are  not  quantitative.  Note  that  “off’  and  “on”  would  work  just  as  well,  but 
low  and  high  are  the  usual  terms. 

There  are  two  methods  for  denoting  a  factor-level  combination  in  a  factorial  design.  The  first  uses 
letters  and  is  probably  the  most  common.  For  example  we  denote  a  factor  level  combination  by  a 
string  of  lower-case  letters  such  as  bed  to  represent  low  level  of  factor  A,  high  level  of  factor  B, 
high  level  of  factor  C  and  high  level  of  D.  The  second  method  uses  numbers  and  generalizes  to 
three-level  and  higher-order  factorials  as  well.  A  factor-level  combination  is  denoted  by  n  binary 
digits,  with  one  digit  giving  the  level  of  each  factor:  a  zero  denotes  a  factor  at  its  low  level,  and  a 
one  denotes  a  factor  at  its  high  level.  Thus  000  is  all  factors  at  low  level,  the  same  as  (1),  and  Oil 
is  factors  B  and  C  at  high  level,  the  same  as  be. 

A  2n  factorial  experiment  was  considered  and  denoted  the  treatments  by  (x1,x2, ...  ,xn ),  where 
xr  =  0  or  1,  r  =  1,2,  ...,n.  and  the  “  true  effect”  of  the  treatment  (x1,x2, ... , xn)  is  denoted  by 
0(x1,  x2, ... ,  xn)  [15].  Then  obtained  the  following  model; 
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0(x1,x2,  -,xn) 

n 

=  F  +  ^  aiFi 
i= 1 

k  n 

+  ^  a^a^F^  +  ^  ai1ai2ai3Fili2i3  +  --  +  a1...akF12^n 

t1,t2  =  l;t1<t2  t1,t2,t3  =  l;t1<t2<t3 

(1) 

Where  is  the  general  mean,  Ft  the  mean  effect  of  the  ith  factor,  Fiii2  the  two  factor  interaction 
between  the  and  i2,  and  so  on,  at  —  1  or  —  1  according  as  Xj  =  1  or  0,  i  —  1,  Let 
y(x1,x2,  —,xn)  be  the  observed  response  corresponding  to  (x1,x2, ...  ,xn).  The  set  of  all  level 
combinations  can  be  represented  by  a  2 n  x  n  Matrix  of  —  l's  and  +l's,  where  +l's  represent  the 
two  levels  of  each  factor,  respectively  see  [16]. 

As  earlier  mentioned  a  convenient  way  of  obtaining  set  of  treatment  combinations  is  to  apply  the 
Yates  algorithm.  Yates  order  state  as  follows:  Write  down  the  observed  treatment  means  in 
standard  order.  With  n  =  3  and  s  =  2,  following  this  order,  the  eight  treatment  combinations  are 
(1),  a,  b,  ab,  c  ac,  be  and  abc,  where  (1)  denotes  the  control,  a  the  application  of  treatment  A 
only,  ab  the  application  of  treatments  A  and  B  and  so  on.  Two-level  full  factorial  and  fractional 
factorial  designs  are  widely  used  in  industrial  experiments  [17]  and  also  in  agricultural 
experiments  [18]  and  [19]  to  assess  the  impact  of  factorial  effects  on  a  process. 

The  estimates  of  the  factorial  effects  is  obtained  on  the  assumptions  that  the  observations  are 
uncorrelated  and  have  equal  variance,  then  in  the  case  of  2n  factorial  designs  provide  independent 
minimum  variance  estimates  of  the  grand  average  and  of  the  2n_1  factorial  effects. 

Consider  a  factorial  design  with  3  two-level  factors.  The  set  of  all  level  combinations  for  the  3 
independent  factors  can  be  written  as: 


/A  B  C  \ 

-1  -1  -1 

-1  -1  1 

-1  1  -1 

D  =  -  111  (2) 

1  -1  -1 

1-1  1 

1  1  -1 J 

\  1  1  i/ 

In  general,  for  n  independent  factors,  the  matrix  D  obtained  in  a  similar  fashion  is  called  the  2" 
full  factorial  design  matrix.  The  set  of  columns  corresponding  to  all  the  main  effects  and 
interactions  is  called  the  2"  full  factorial  model  matrix,  denoted  by  [16]. 
where  X  —  2”  x  (2"  - 1)  matrix  defined  by 

Hn  =  {A, B,  AB,C,  AC, BC,ABC,...,ABCD...k}2„  and  n  —  1,2, .../c  with  entries  -1  and  1  and 

columns  arranged  in  Yates  order,  where  A  =  (-1,...,— l,l,...,l) , 

B  =  (-1,.. .,-1,1,. ..,1,-1,. ..,-1,1,.. .,l)  ,..., and  n  =  (— 1,1,— 1,1,..., — l)  stand  for  its  n  independent 

columns,  and  the  other  columns  are  expressed  as  products  of  the  n  independent  columns. 

The  corresponding  model  matrix  Xfor  the  2 3  design  is  given  by 
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X  = 


(A  B  C  AB  AC  BC  ABC\ 
-1-1-1  1  1  1-1 

-1-1  1  1  -1  -1  1 
-1  1-1-1  1-1  1 

-1  1  1-1-1  1  -1 
1  -1  -1  -1-11  1 
1-1  1-1  1-1  -1 


(3) 


This  representation  of  the  factor  level  combinations  is  convenient  since  the  columns  of  X  denote 
the  linear  contrasts  that  estimate  the  main  effects  and  interactions  in  a  normal  linear  regression 
XY 

model  by  - ,  where  Y  is  the  vector  of  observations  corresponding  to  the  factor  level  settings  of 


each  row  of  D. 


1.2  Lexicographical  Order 

Lexicographical  method  is  similar  to  dictionary  (alphabetical)  order  [20].  If  Word  A  is  shorter 
than  Word  B.  and  every  letter  of  Word  A  occurs  in  the  same  place  in  Word  B.  Word  A  comes 
before  Word  B  (’’compute”  and  ’’computer”)  If  the  first  letter  that  differs  in  Word  A  comes  before 
the  corresponding  letter  in  Word  B,  then  Word  A  comes  before  Word  B  (’’math”  and  ’’matter”). 

2.  Methodology 

2.1  Some  Proposed  Methods  for  Generating  Factorial  Effects  and  their  Treatment 
Combinations 

The  classical  method  of  generating  factorial  effects  is  inefficient  of  time  and  space.  Hence  we 
consider  a  new  method  that  leads  to  a  more  formal  and  simpler  in  order  to  obtain  main  effects, 
interactions  and  treatment  combinations. 

Some  alternative  methods  for  generating  factorial  effects  and  their  treatment  combinations  are 
proposed.  These  include  using  the  general  concept  of  Pascal  triangle  and  lexicographical  order. 

2.1.1  A  New  Method  of  Generating  2”  Factorial  Effects  from  Pascal  Triangle 

Under  this  subheading  a  proposition  was  made  to  obtain  both  the  number  of  factorial  effects  and 
their  treatment  combinations  as  follows: 

Proposition  2.1  The  generation  of  orthogonal  contrasts  in  2n  factorial  design  follows  Pascal 
triangle. 

Proof. 

In  2"  where  n  ~  be  the  real  numbers.  Then  prove  by  induction. 

For  ^  ~  -  the  result  is  trivially  true. 

Now  it  is  also  true  for  h  =  k  since  we  have: 
n  —  main  effects, 

n(n  —  1)  n(n— l)(n  — 2) 

-  two  —  factor  interaction, -  three 

2!  1  2x3 

factor  interaction  effects, ..., 71(71  — —  ^  ^  k+1'>  ~  (n  —  1)  factor  interactions  and  finally 
In  —  factor  interaction. 
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Finally  for  h  =  k  +  1 ,  we  obtain  n  -  main  effects, 


n(n- 1)  n(n-l)(n-2) 

2!  two  _  factor  interactions,  2x3 


three  -  factor  interactions, 


n(n-\)(n-2). .  .(u-(k+l)+l)  n(n-\)(n-2). .. tyi-k ) 


(k+1)! 


(k+1)! 


Let  N  =  k  +  l 
Jfc  =  JV- 1 

n(n-l)(n-2)...(n-k)  n(n-\)(n-2)...(n-(N-l)) 


n(n-l)(n-2)...(n-N+l) 


which  is  also  true  for  h  =  k  +  1 . 


(k+1) !  N !  N ! 

2.2  Lexicographical  Method  of  Generating  Treatment  Combinations 

Without  loss  of  generality,  let  us  assume  that  the  v  =  sk  treatment  combinations  for  sk  where 
s  =  slf  s2, ... ,  (s  >2)  and  k  >2  this  can  be  generated  lexicographically  by  the  following  rule: 
Rule  I :  This  rule  is  applied  to  S  levels  at  2  factors  ( i.e  S=2,3,...  and  k=2) 

1.  Initialize  the  first  treatment  combination  of  all  factors  at  lower  level. 

2.  Increase  the  level  of  2nd  factor  holding  other  factor  at  lower  level  up  to  s~2 

3.  Repeat  step  2  but  increasing  the  level  of  first  factor  up  to 


Rule  II:  This  rule  is  applied  to  S  levels  at  3  factors  ( i.e  S-2,3,...  and  k=3) 

1.  Initialize  the  first  treatment  combination  of  all  factors  at  lower  level. 

2.  Increase  the  level  of  2nd  factor  holding  other  factors  at  low  level  up  to  s~2 

3.  Repeat  step  2  but  increasing  the  level  of  first  factor  up  to  s[. 

4.  Taking  into  consideration  of  the  previous  steps  allow  the  level  of  the  third  factor  to  vary  up 
to  ^ 

Were  si  =  sx  —  1,  si  =  s2  —  1  and  s~3  —  s3  —  1. 


3.  Results  and  Discussion 

3.1  Application  of  Pascal  Triangle  in  Generating  Number  of  Factorial  Effects 

The  classical  method  of  generating  factorial  effects  lead  to  spending  of  time  and  space.  Hence  the 
need  to  consider  a  new  method  that  leads  to  a  more  formal  and  simpler  in  order  to  obtain  the 
number  of  main  effects  and  interactions.  From  the  Proposition  2.1  for  generating  factorial  effects 
from  Pascal  triangle  as  earlier  stated.  We  obtain  the  following  in  the  case  of 
2"  where  n=0, 1,2,3, ...,k 


The  Pascal  triangle  for  n  is  written  as  follows: 


1 


1 


1 

1 

5 


1 


4 


1 


2 

3 

6 


10 


_ ^  n  —  0 

1  - ►  n  =  1 

1 - ►  n  =  2 

3  1  - >n  =  3 

4  1 - ►  ?!  -  4 

10  5  1 - ►  n  =  5 


1 

2 

4 

8 

16 

32 


1 


k 


k(fc-l)  fc(fc-l)(fc-2) 

2!  3! 


k  1— ►  n  =  k  2k 
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For  n  =  2 

1  2  1 

1  1  1 

1  —  identity  2  —  me  1  —  2fi 

This  shows  that  in  22  factorial  design  the  triangle  generates  1 —mean  (/)  the  identity, 
2— main  effects  (me),  and  1—2— factor  interactions  (fi). 

For  n  =  3 

1  3  3  1 


1  —  identity  3  —  me  3  —  2  fi  1  —  3  fi 

For  23  factorial  design  the  triangle  will  generates  1-mean  (/)  the  identity,  3- main  effects  (me) 
,  3— 2— factor  interactions  (fi),  and  1—3— factor  interactions  (fi). 


For  n  —  4 


b  4  1 

1  1  1 


1  —  identity  4  —  me  6  —  2 fi  4  —  3 fi  1  —  4 fi 

For  24  factorial  design  the  triangle  will  generates  1 —mean  (/)  the  identity,  4— main  effects  (me) , 
6— 2— factor  interactions  (fi),  6— 2— factor  interactions  (fi),and 
1—4  —  factor  interactions  (fi). 

For  n  =  5 

5  10  10  51 


'f  t  ▼ 


1  —  identity  5  —  me  10  —  2  fi  10  —  3  fi  5  —  4  fi  1  —  5  fi 

For  25  factorial  design  the  triangle  will  generate  1 —mean  (/)  the  identity,  5— main  effects  (me), 
10— 2— factor  interactions  (fi),  10—3 —factor  interactions  (fi), 

5— 4— factor  interactions  (fi), and  1—5— factor  interactions  (fi). 

For  n  —  k 


K 

1 


fe(fc-l) 


fc(fc-l)(fc-2) 


2! 


3! 


1  —  identity  k  —  me 


fc(fc-i) 

2! 


-2  fi 


fe(fe— l)(fc— 2) 
3! 


—  3fi  k  —  (k  —  1  )fi 


t 

1  —  kfi 


For  2  factorial  design  the  triangle  will  generates  1  —mean  (/)  the  identity,  k—main  effects  (fi) , 

k(k-l)  .  .  .  k(k-l)(k-2)  0  .  . 

— — - 2- factor  interactions  (h),  - — - 3- factor  interactions  (h),  , 

k—(k—l)— factor  interactions  (fi),  and  1  -k— factor  interactions  (fi). 
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3.2  Application  of  Lexicographical  Method  in  Generating  the  Treatment  Combinations 


The  above  rules  I  and  II  to  are  being  applied  to  obtain  the  following  table  for  cases  of 
(s  =  2,fc  =  2),(s  =  2,fc  =  3),(>  =  3,fc  =  3)and(s  =  3,fc=4). 

Table  1.  The  List  of  SK  Treatment  Combinations  with  S  prime  number  and  K  number  of  factors 


Where  5,  =  sl  - 1 ,  s2  =  ,s2  - 1 ,  s3=s3-land  s4=s4—l 

Table  1  display  the  treatment  combinations  of  the  SK  with  the  algorithm  of  how  they  are 
generated  using  K  —  2,3A  and  S=2,3-  The  rigorous  mathematical  computation  involved  in 
kronecker  product  formulation  of  [2]  in  generating  factorial  effect  and  treatment  combinations 
make  it  hard  and  cumbersome.  The  proposed  method  is  a  convenient  one  in  achieving  the  same 
target.  While  [2]  restrict  the  discussion  on  symmetric  2-level  factorial  designs  this  paper  consider 
both  symmetric  and  asymmetric  s-level  factorial  design  as  can  be  seen  in  the  application  of  the 
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algorithm  for  lexicographical  order.  Also  the  proposition  under  Pascal  triangle  can  be  used  to 
calculate  number  of  combinations  of  factorial  effects  and  treatment  combinations  using  the 
recursion  formula. 

4.  Conclusion 

The  classical  method  of  generating  treatment  combinations  is  only  applicable  to  Yates’  forward 
algorithm  in  estimating  the  effects  of  full  factorial  and  the  aliased  subsets  of  fractional  design.  The 
Yates’  order  which  is  difficult  to  implement  in  fractional  factorial  design  especially  in  the  case  of 
n  >  3 .  This  is  mostly  common  in  generating  the  blocks  elements  of  principal  block. 

The  existing  method  will  pave  way  for  estimating  the  factorial  effects  for  not  only  the  full  and 
fractional  factorial  designs  but  for  designs  with  confounding,  partial  confounding  and  split  plot. 
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